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LOCAL AFTERBODY HEAT TRANSFER TO A BLUNT 


TWO-DIMENSIONAL CONFIGURATION AT MACH 8* 

By Dennis M. Bushnell 
Langley Research Center 

SUMMARY 

Local heat-transfer rates and pressures were measured on a blunt two-dimensional 
configuration the afterbody of which was subjected to an extensive region of separated 
flow. The investigation was conducted at a free- stream Mach number of 8 and at free- 
stream Reynolds numbers based on the model face width of 2 inches (5.08 cm) from 
1.27 X 10 5 to 7.55 X 10 5 . 

The predictions from two theories were compared with the heat-transfer data 
obtained in the afterbody region. One of these theories was a constant-pressure integral 
method utilizing profiles from similar solutions for enthalpy and velocity distributions. 
The other theory was a local flat-plate method in which the velocity and enthalpy condi- 
tions obtained from the integral method in the reverse-flow region were used as effective 
free-stream quantities. Both theories were in reasonable agreement with the general 
level of the data. The integral method predicted the correct trend in the vicinity of the 
separation point, whereas the flat-plate method gave the correct trend over the aft por- 
tion of the afterbody. 


INTRODUCTION 

The current use of atmospheric braking during the reentry of a space vehicle has 
necessitated the development of methods for predicting the severe heating which occurs 
during deceleration. For the forward portion of the vehicle, where the flow is attached 
and the heat transfer is the highest, theoretical predictions have been successful. How- 
ever, for the aft portion of the vehicle, which is typically a region of separated flow, a 
satisfactory method of predicting local heating has not yet been developed. The flow 
separation which occurs over the aft portion of a vehicle is usually termed "free" sepa- 
ration, as opposed to the "forced" separation which is generally caused by deflected flaps 
or shock impingement. 

*The information presented herein was offered as a thesis in partial fulfillment 
of the requirements for the degree of Master of Science in Mechanical Engineering, 
University of Virginia, Charlottesville, Virginia, June 1967. 



Several approximate methods are currently available for predicting the local heat 
transfer in a "forced" laminar separated flow which reattaches on a surface. All of 
these methods employ some integral form of the boundary -layer equations. The first 
method appeared in reference 1, where the assumptions of polynomial profiles for the 
velocity and static enthalpy distributions, constant pressure, and zero boundary -layer 
thickness at the separation point were made. Primarily because of the latter assump- 
tion, the method is not applicable in the vicinity of the separation point. 

The second method (ref. 2) is an approach similar to earlier work in references 3 
and 4, except for the addition of the integral energy equation. This method employs the 
reverse-flow velocity and enthalpy profiles from the similar solutions in reference 5. 

It was successfully used in predicting local heat transfer in separated flows caused by 
deflected flaps and shock-wave boundary- layer interactions. 

A third approach to the problem of predicting local heat transfer in a region of 
forced separation was presented in reference 6, where the Dorodnitsyn method of integral 
relations was employed. Although the method could be used to predict local heat trans- 
fer, no calculations or comparisons with experimental data were presented. 

Because the second method (ref. 2) has thus far been successful in predicting local 
heat transfer for forced separation, it may be expected to provide good predictions for 
the free- separation problem. The purpose of the present investigation is to compare 
with experimental afterbody data the results of an approach similar to that of reference 2 
as applied to afterbody flows, but simplified by the assumption of constant pressure (see 
appendix A for details of method). The data used for comparison were obtained at a 
Mach number of 8 on a two-dimensional blunt configuration which had an extensive region 
of separated flow at essentially constant pressure. 


a 

b 

c l ,c 2 ,c 3^l 
c 4> c 5’ c 6 J 


c p 

D 


SYMBOLS 
Y —0 

velocity profile parameter, - u ~ - 
enthalpy profile parameter, 

coefficients in polynomial curve fit of integral parameter 
(tables I and II) 

specific heat at constant pressure 

model face width, 2 inches (5.08 centimeters) 
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dj,d 2 ,d 3 > j coefficients in polynomial curve fit of integral parameter 
(table III) 
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function related to stream function (eq. (23)) 

„2 

stagnation enthalpy, h + 
static enthalpy 
integral parameter, 
integral parameter, Q*j 6* 

thermal conductivity 

distance along afterbody of model (see fig. 1(a)) 
Mach number 
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integral parameter, 5^ 
pressure 
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heat-transfer rate 
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free- stream Reynolds number based on model face width 
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skin thickness 
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temperature 
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integral parameter, 
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6*(tw - 1) 
velocities in x- and y- directions 

velocities in transformed coordinate system (eq. (8)) 
longitudinal and normal coordinates 

transformed longitudinal and normal coordinates (eq. (5)) 
distance defined in equation (15) 
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similarity variable (eq. (24)) 
angle defined in figure 1(a) 
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parameter defined in equation (6) 


dynamic viscosity 


kinematic viscosity 
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p density 

t time 

<p = -&±- 

PwM w 

ip stream function (eq. (7)) 

Subscripts: 

e local flow quantity external to boundary layer 

o separation point 

rev conditions at maximum velocity in reverse-flow region 

s local stagnation conditions 

t stagnation conditions behind a normal shock at free-stream conditions 

w wall conditions 

00 free-stream conditions 

* value at dividing streamline 

Single, double, and triple primes denote first, second, and third derivatives with 
respect to 77 . ^Note that, in addition, f* = ~ 

APPARATUS 

Tunnel 

The model tests were conducted in the Langley Mach 8 variable-density hypersonic 
tunnel. This tunnel has an axially symmetric nozzle with contoured walls. The test 
medium is air. The average test- section Mach number variation with stagnation pressure 
is given in reference 7. For the present investigation the stagnation pressure was varied 
from 130 to 1055 psia (0.896 to 7.27 MN/m^). The corresponding Mach number range is 
from 7.74 to 7.95. The range of free-stream Reynolds numbers based on the 2-inch 
(5.08-cm) model face width was from 1.27 x 10^ to 7.55 X 10^. Nominal stagnation tem- 
perature was 900° F (755° K). Further details concerning the facility are available in 
reference 8. 
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Model 


General descripti on.- A sketch of the model is shown in figure 1(a). The configu- 
ration is basically a section of a cylinder which was mounted normal to the free- stream 
flow direction with the blunt face forward at zero angle of attack. The supporting struts 
attached to the ends of the model did not extend above the model center line so that 
schlieren photographs of the afterbody flow could be taken. The model was constructed 
with a thin skin silver- soldered to internal supports. Type 347 stainless steel was used, 
and the model had a large ratio of span to face width (4.5) in order to minimize end 
effects. The nominal skin thickness was 0.030 inch (0.076 cm). Photographs showing 
front and oblique views of the model are presented in figure 1(b). 

Instrum entation . - Figure 1(a) also shows the location of the pressure orifices and 
thermocouples. Pressure orifices were installed on the model at the indicated positions 
with the use of monel tubing having an inside diameter of 0.040 inch (0.102 cm); the 
tubing was mounted flush with, and normal to, the local surface. The afterbody of the 
model was also instrumented with thermocouples made from No. 30 chromel-alumel 
wire 0.010 inch (0.025 cm) in diameter. The thermocouples were installed in the posi- 
tions shown by drilling two 0.011-inch-diameter (0.0279-cm-diameter) holes 0.020 inch 
(0.0508 cm) apart, inserting the individual thermocouple wires into these holes, silver 
soldering the wires in place, and fairing the surface flush. 

Equipment 

The temperature-time history of the model was recorded on magnetic tape with an 
analog-to-digital data recording system at a rate of 40 samples per second. Each pres- 
sure orifice on the face was connected to two strain-gage pressure transducers - one 
with a range from 0 to 1 psia (0 to 6.89 X 10^ N/m2) and the other with a range from 0 
to 10 psia (0 to 6.89 x 10^ N/m^). The orifices on the afterbody were connected to 
thermocouple vacuum gages with a nominal range from 0.01 psia (6.89 x 10 N/m2) to 
0.3 psia (2.07 X 10^ N/m2). The output from these pressure-measuring devices was 
also recorded by the analog-to-digital data recorder at 40 samples per second. Tunnel 
stagnation pressure was measured with a Bourdon tube gage. 

TEST PROCEDURES AND DATA REDUCTION 
Heat- Transfer Data 

The heat-transfer data were obtained by the transient heating technique. For this 
technique, steady flow was established in the test section with the model retracted from 
the flow; the model was then injected into the test section, the total injection time being 
approximately 0.2 second. Data were taken for approximately 8 seconds, after which 
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the model was removed from the airflow and brought to isothermal conditions at approxi- 
mately room temperature (90° F (305° K)) in preparation for another test. 

The temperature data were reduced to heating rates by use of the calorimeter 
equation 

q = PCpS (1) 

wherein both conduction and radiation heat transfer are neglected. The radiation heat 
transfer was measured by inserting the model into the tunnel test section with no airflow 
and monitoring the skin- temperature variation with time. From such tests the radiation 
heat transfer to the model from the test- section walls was determined to be approxi- 
mately 0.5 percent of the lowest aerodynamic heating rate measured and was therefore 
negligible. 

For test times of less than 1 second, when the model was essentially isothermal, 
the maximum heat transfer by conduction along the skin at the thermocouple locations 
was calculated to be approximately 8 percent of the lowest aerodynamic heating rate 
measured. As this conduction error was generally smaller than other experimental 
errors, a conduction correction was not applied during the data reduction and equation (1) 
was used to calculate the aerodynamic heating. The local skin thickness was measured 
for each thermocouple position. The values of p and Cp used were 500.5 lbm/ft^ 

(8 x 10^ kg/m^) and 0.12 Btu/lbm-°R (502 joules/kg-°K), respectively (ref. 9). In addi- 
tion to conduction effects there are possible errors in measured temperatures, skin 
thickness, and the values of density and specific heat of the skin; hence, the final values 
of heating rate are considered accurate to within 15 percent. 

Pressure Data 

The pressure data were obtained during the same runs as the heat-transfer data. 

The run time of 8 seconds was sufficient to allow the gages to reach a steady-state con- 
dition. As the measured wall temperature on the afterbody did not change more than 
20° R (11.1° K) during this interval, the pressure and heat- transfer data were obtained 
at essentially the same value of wall temperature. The pressure sensors were calibrated 
before the tests and checked during the testing period. The transducers that were used 
to measure the pressure over the face of the model are accurate to within 1 percent of 
the full-scale reading. Therefore, the maximum error in the face-pressure data is 
5 percent. The thermocouple vacuum gages used for the afterbody pressure readings 
are accurate to within 0.5 percent of the full-scale reading. However, because of the 
low pressures on the afterbody and the difficulty of calibration in this range, the final 
data are accurate only to within 10 percent. Schlieren photographs were obtained during 
the same tests as the heat-transfer and pressure data. 
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RESULTS AND DISCUSSION 


Flow Field 

Shown in figure 2 are schlieren photographs of the flow field over the model for 
several values of free-stream Reynolds number. The dark area in the vicinity of the 
model corner indicates an expansion region. The light area immediately beneath this 
expansion region indicates a compression region, and therefore the upstream flow was 
apparently overexpanded. (The expansion and compression regions under discussion 
are labeled in fig. 3(a).) To ascertain if the overexpansion was large enough to cause 
the flow to attach to the afterbody, this section of the model was coated with oil and then 
injected into the airstream. For a small distance just downstream of the corner, the 
oil was removed by the airstream (see fig. 3(b)), an indication that the flow does remain 
attached for this distance along the afterbody. The fact that the oil was undisturbed on 
the remainder of the afterbody indicated the presence of a low -shear region of separated 
flow. The attachment of the flow for some small distance downstream of the corner is 
in agreement with results reported in reference 10 for flow separation from a flat-based 
cone- cylinder configuration. 

The dark band which appears underneath the compression region on the photographs 
in figures 2 and 3(a) at an angle to the afterbody represents a region of decreasing den- 
sity as the body is approached (inasmuch as the schlieren method gives density gradient). 
Since the pressure is nearly constant in this region, this decreasing density corresponds 
to an increase in temperature as the flow velocity is reduced. The largest temperature 
gradient generally occurs in the outer portion of the viscous mixing region, which grows 
after flow separation occurs; hence, the dark band corresponds to this portion of the 
flow region. The inclination of this dark band with respect to the free-stream flow direc- 
tion becomes greater with increasing Reynolds number, an indication that the ratio of the 
afterbody pressure to free-stream pressure should decrease with increasing Reynolds 
number. For the purpose of calculating afterbody pressure, the slope of the dark band 
under discussion was assumed to be equal to the local flow angle at the edge of the vis- 
cous region. The location at which this slope was measured is shown in figure 3(a). On 
some of the photographs in figure 2 (ones for the higher values of R^d)* the neck of 
the wake and the beginning of the far wake are visible. 

Pressure Distributions 

Face .- Figure 4 shows the measured pressure distribution over the model face 
for the range of free-stream Reynolds numbers of this investigation. As expected, the 
distribution over the face is essentially independent of Reynolds number. Also shown 
in figure 4 is a modified Newtonian distribution. The surface pressure falls below the 
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Newtonian value as 6 increases, as previously noted for sections of spheres (see 

ref. 11). In the heat-transfer theory summarized in appendix A, the quantity ^ P e u e 

is required, with the integral evaluated from the stagnation point to the separation point. 
To evaluate this quantity over the model face, the dashed line shown faired through the 
data in figure 4 was used. This curve was faired into the sonic pressure ratio at the 
corner (0.528). 

Afterbody .- Figures 5(a) to 5(c) illustrate typical measured afterbody pressure 
distributions. The figures also show the level obtained by assuming an isentropic 
expansion from the stagnation point to the afterbody flow angle measured in accordance 
with the method described in the flow-field discussion. For this calculation, the flow 
angle for the Mach 1 condition was assumed to be equal to the flow angle behind an 
oblique shock which just gives a Mach number of 1 downstream of the shock at the free- 
stream conditions. This angle is approximately 43.6° and was obtained from the tables 
of reference 12. (The Mach 1 inclination near the surface of the model would be greater, 
as indicated by the data in figure 4.) The flow was then assumed to undergo a Prandtl- 
Meyer expansion from 43.6° to the final flow angle measured on the afterbody from 
schlieren photographs like those of figure 2. As discussed previously, this flow angle 
was taken as the slope of the dark band at the point indicated in figure 3(a). In order to 
calculate the afterbody pressure by the method described, it was necessary to assume 
that the recompression of the afterbody flow discussed in the flow-field section was 
isentropic. 

The other pressure level shown in figure 5 was calculated by using the method of 
reference 13. The value of u*/u e used in the calculation was 0.587 which is the value 
given by the similar solution in reference 14 for an external stream mixing with an 
infinite quiescent reservoir. In reference 3 it is shown that this value will usually be 
achieved only far downstream in an actual laminar mixing region and that for low values 
of Reynolds number this value may never be reached in the near wake. Therefore the 
value of u*/u e of 0.587 and the method of reference 13 should give a value of afterbody 
pressure which would be approached asymptotically as the Reynolds number is increased 
(with the assumption that the afterbody flow remained laminar). This phenomenon is 
discussed in more detail in reference 15. 

In figure 5 the data indicate that the afterbody pressure is essentially constant. 
Also, the general level of the experimental data is in reasonable agreement with the 
level calculated by using the measured afterbody flow angle. In figure 5(a), which is for 
a low Reynolds number, the level calculated by the method of reference 13 is consider- 
ably below that of the data; however, as expected, the data at higher values of Roo^d are 
in closer agreement with this calculated level (see figs. 5(b) and 5(c)). 
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To show the variation of afterbody pressure with Reynolds number, the pressure 
at three values of l / D has been plotted as a function of Reynolds number in figure 6. 
From this figure it is evident that as the Reynolds number is increased up to about 
4 x 10 5, the afterbody pressure ratio tends to approach the value predicted by the method 
of reference 13. For R,*, p >4x 10^ the pressure ratio drops below the value pre- 
dicted by this method. This latter result may be caused by transition to turbulent flow 
in the near wake. The near -wake transition studies for axisymmetric blunt bodies 
reported in reference 16 (wherein a model cross section somewhat similar to the one 
used in the present investigation was employed) indicate that 4 x 10 5 would be a reason- 
able Reynolds number at which to expect transition in the near wake. 

Heat- Transfer Distributions 

Variation along aft erbody.- The measured afterbody heat-transfer distributions 
are presented in figure 7, each part of which represents a different value of test Reynolds 
number. Also shown in the figure is the distribution calculated by using the method 
given in appendix A and the measured pressure level. The method (appendix A) is essen- 
tially that of reference 2 except that it is simplified for application to the afterbody 
problem by the assumption of constant pressure. The integral constant-pressure 
boundary -layer equations are solved by using the reverse-flow profiles of reference 5, 
which are uncoupled from the pressure- gradient parameter. The pressure level used 
in the calculation is noted in each figure. 

The value of % for the actual data varied from 0.386 to 0.44. Inasmuch as the 
theoretical calculations were carried out for a value of tw (wall- to total -temperature 
ratio T w /T s e ) of 0.35 to obtain the distribution shown in figure 7, it was assumed that 
the variation of %w/{^w)x-0 X (fig. 8) was independent of t w over the interval 

from 0.35 to 0.44. The value of (?^)~ q f° r the t w value pertinent to the test under 

consideration was used to obtain as a function of X. The pertinent value of t w 
is given in each part of figure 7. 

As shown in figure 7(a), the theory of appendix A indicates a large decrease in 
heat-transfer rate near the separation point, followed by a more gradual decrease along 
the rest of the afterbody. The data tend to agree with the prediction of the initial 
decrease, but as l increases the data show a trend of increased heating. 

The distribution of the data appears reasonable if the physical flow picture is con- 
sidered. The flow adjacent to the afterbody surface is a region of reverse flow. This 
flow might be expected to cause a sublayer (roughly equivalent to a flat-plate boundary 
layer) to form along the body. This equivalent flat-plate boundary layer would start at 
the rear of the body and grow in the negative l direction. Therefore, the heat-transfer 
rate might be expected to decrease as l decreases. The theory of appendix A 
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apparently fails to give the proper heat-transfer-rate distribution because this boundary 
condition at the rear of the model cannot be included in the theory. 

This analogy between the sublayer near the afterbody surface and a flat-plate 
boundary layer can be utilized to calculate the afterbody heat-transfer rate by an approxi- 
mate method that may be more realistic away from the separation point than is the 
method of appendix A. Basically, the inner (reverse-flow) portion of the velocity and 
enthalpy profiles is assumed to constitute an effective free-stream flow, which causes 
a flat-plate boundary layer to begin growing in the negative l -direction at the rear of the 
model afterbody. The nominal values of f' and £ that characterize this flow were 
taken from the profiles obtained in the vicinity of the rear of the model by the method of 
appendix A. These values were f^. ev = -0.05 and ?rev = 0*5 and are characteristic 
of the fairly flat maximum-velocity portion of the reverse-flow region. The value of 
as computed from the results of reference 17 for /3 = 0 (flat-plate) attached flow 
is 0.305. With this information and the local measured pressure level, a flat-plate heat- 
transfer- rate distribution was calculated by the method of reference 17 and is shown in 
figure 7 as the dashed lines. (The rear of the model was taken as the origin of the flat- 
plate flow.) In all parts of figure 7, flat-plate theory seems to predict the distribution 
away from the separation point, except in the immediate vicinity of the rear of the model. 

In summary, the theory of appendix A correctly predicts the approximate level and 
the measured trend in heating near the separation point but not over the rest of the after- 
body. The flat-plate theory predicts essentially the correct type of distribution over the 
rest of the afterbody, but the predicted level of heating does not increase as much with 
Reynolds number as the data. 

After the present investigation had been completed, an approach to the prediction 
of afterbody heat transfer that was similar to the flat-plate method used herein appeared 
in reference 18. The approach was to solve the "free- shear- layer" problem and then 
apply the blunt-body attached-flow heat-transfer expression to the afterbody with the 
average separated- region quantities as effective external conditions. Base heating to a 
flat-based ogive- cylinder was successfully predicted by using this method. These 
results, together with those of the present investigation, indicate that the use of attached- 
flow expressions with local recirculating or reverse-flow quantities as effective external 
conditions appears to be a valid approach to the prediction of afterbody heat transfer in 
regions of separated flow. 

Variation with Reynolds number .- In order to show a comparison between the 
experimental variation of heat transfer with Reynolds number and the variation predicted 
by the flat-plate theory, the heat-transfer rate at l/T> = 0.47 is plotted as a function of 
the local unit Reynolds number (p e u e/ /ju e ) in figure 9. It is apparent that the increase in 
heat- transfer rate with unit Reynolds number for p e u e /p e > 1.5 x 10 ^ per foot 
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(4.92 X 105 per meter) is greater for the data than for the theory. This increase with 
unit Reynolds number is probably associated with transition to turbulent flow in the near 
wake, as the value of unit Reynolds number where the increase occurs corresponds to 
Roo d = 4x 105, which is the "transition Reynolds number" indicated by the pressure 
data of figure 6. (See discussion of fig. 6.) 


Engineering estimate s of afterbody heat transfer.- By using the flat-plate method, 
it may be possible to obtain approximate engineering estimates for heating in a region of 
separated afterbody flow without solving the equations of appendix A to obtain values for 
*rev anc * ?rev indicate a method by which these estimates might be obtained, fig- 
ure 10 shows the quantities f^ev anc * Crev i rom reverse-flow similar solutions for 
Np r _i and <p = 1 plotted as functions of the velocity and enthalpy profile parameters 
at several values of t^. 


Figure 10(a), in which the values of C rev have been normalized by using the 
appropriate value of shows that the normalization procedure partially collapses the 
curves for different values of t w , but overemphasizes the variation of ?rev with b 
for a given tw Nevertheless, the fact that the curves are fairly flat and tend to group 

about a value of — — of 0.2 indicates that this value would provide a reasonable 

estimate for the reverse-flow enthalpy parameter over a wide range of values of b and 
tw; therefore, a solution for b by the method of appendix A may not be critical to 
obtain a reasonable estimate for £ rev . 


The values of f^ ev shown in figure 10(b) exhibit a large variation with the param- 
eter a. Therefore, a local estimate for a is needed in order to apply the flat-plate 
method. For large values of a (thick separated regions), an upper heating limit could 
be provided by using values for f^ ev ranging from -0.09 to -0.12 depending on the 
value of tw- 


For smaller values of a (less than 0.4) and in the absence of a solution to the 
equations of appendix A, an estimate for a can be made by extrapolating the initial 
slope of a with X. This initial slope can be obtained by solving equations (16) and (19) 
of appendix A for da/dX. Note that dH/dX = (dH/ da) (da/ dX) and, for example, at 
X = 0 (as the initial a slope is being ascertained) P = 0, 6? = 1, and a = 0. The 

result is 



R 


J dH + d J 
H da da/ a=0 


The values of the integral parameters can be obtained from the coefficients in 
table I for t w = 0.35 and from reference 4 for V = 1 and t w = 0.2. The extrapola- 
tion of the initial slope gave a value for a near the rear of the model about 50 percent 
lower than the actual value computed by the method of appendix A; the heating would 
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therefore be underestimated by about 60 percent, as compared with the results obtained 
when using the method of appendix A. Thus, this extrapolation method is fairly crude 
and should be used only if rough estimates are desired. 

Another possible method for estimating a value of a with which to enter fig- 
ure 10(b) is to use an estimate of the physical distance to the dividing streamline. The 
procedure would be to obtain an afterbody flow angle (assumed constant) from correla- 
tions or knowledge of the afterbody pressure. The distance from the body surface to the 
local flow line drawn from the separation point parallel to the afterbody flow direction 
would then provide a rough estimate of the physical distance to the dividing streamline. 
From equation (24) 


C v * p e , 
7*i = 


Jo 




2 J 0 p w p w u e * 


i 


When the information on the right-hand side of this equation is known, a value of 
V* p e 


_ d 77 with which to enter the reverse-flow velocity profiles can be obtained. 

0 P 

The integral over x ' should be taken from the stagnation point. The value of a asso- 
rt* P e 

ciated with the profile having the known value of \ — dri is then an approximation 

J 0 P 

of the local value of a. 

In general, for a perfect gas 

H 


r 

J o 




u; 


rv* rV* t 9 

i C d„ - r (f ) dn 


s,e/J 


C 7 /* r 1 '* 2 

The integrals \ £ dri and \ (?) dr] are given in figure 11 as functions of a. An 

J 0 J 0 

approximate 20 percent overestimation of the local a near the rear of the model 
results when this iterative method is applied to the present data by using a value of y* 
(distance to the dark band) obtained from the schlieren photographs of figure 2. Com- 
pared with results obtained when using the method of appendix A, the heating would be 
overestimated by about 15 percent. Therefore, this method appears to give better 
results than the previously discussed procedure of extrapolating the initial slope. 


CONCLUSIONS 


Local heat-transfer rates and pressures have been measured over the afterbody of 
a blunt two-dimensional configuration at a Mach number of 8. The afterbody of the model 
was subjected to an extensive region of separated flow. The model was tested at a 
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nominal value of wall- to total-temperature ratio of 0.4 and at free-stream Reynolds 
numbers based on face width from 1.27 x 10^ to 7.55 x 10^. The heat-transfer data were 
compared with the results of a constant-pressure integral method employing reverse- 
flow profiles and with the results of a flat-plate method. The following conclusions can 
be made: 

1. The flow remains attached over the afterbody downstream of the corner of the 
model for a small distance before separation occurs. 

2. The pressure level is essentially constant over the afterbody and decreases with 
increasing Reynolds number. 

3. The measured heat-transfer rate decreases just downstream of separation and 
then increases in the downstream direction as the rear of the model is approached. The 
integral method correctly predicts the initial large decrease in heating, but also predicts 
a gradual decrease in heating toward the rear of the model, a trend opposite to that indi- 
cated by the data. 

4. Away from the separation point, the flat-plate method predicts the correct type 
of heat-transfer-rate distribution for the lower Reynolds number data. 

5. The deviation from the theoretical laminar predictions of both the heat-transfer 
and pressure data at the higher Reynolds numbers probably results from afterbody-flow 
transition beginning at a free-stream Reynolds number based on face width of about 

4 X 10 5 . 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 16, 1967, 

129-01-08-38-23. 
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APPENDIX A 


DETERMINATION OF LOCAL LAMINAR HEAT- TRANSFER DISTRIBUTION 
IN A CONSTANT-PRESSURE SEPARATED FLOW 


As in the analysis of reference 2, the starting point of the method is the boundary- 
layer equations for a compressible perfect fluid. For the case of no streamwise pres- 
sure gradient, these equations are as follows: 


Continuity: 

Momentum: 

Energy: 


i(pu)+i(pv)=° 


S + w = IK 1 * w) 


The introduction of a modified Lees transformation 




Pe u e 


-\ 


X = \ X dx 

J 0 Poo U oo 


Y = ^r y ^dy 

Uoo Jq P oO 


where 


M oo Too 


and definition of a stream function 4/ 


9 4 s £U_ 

9y Poo 

M = .fix. 


3x 


Poo, 


so that 


u = -^ 

9Y 

y = - Ml 

9X 


( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
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transform equations (2) to (4) to the following analogs of the incompressible form of the 
boundary- layer equations: 

Continuity: 


Momentum: 


Energy: 


U 


£U + 9V = o 
ax ay 

(9) 

9 u + v 9U = a 2 u 

ax ay 00 9y 2 

(10) 

as , „ as a 2 s 

ax v ay ay 2 

(ID 


These equations apply for a perfect gas with constant specific heats and a constant 
Prandtl number of 1 where 


S =-^-- 1=?- 1 

^Sje 


Integration of equation (10) across the boundary layer and the use of equation (9) give 


d0 j _ V arj /9U \ 
dX jj 2 \9Y/ W 


(12) 


Multiplication and division by 6^ give 


'K*) 


where 


and 


dX Ugfij* 


- 0i 


(13) 


P = 5? 


<U/U e ) 


9Y 


Jw 


From the definitions 


and 


er 

= 5. 


Wo 1 


~ X. 

X = 


( 5 i)o] U e 


(14) 


(15) 
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the final form of the momentum equation can be written as 


H 


T* 


(6*f® = P 

' w dX 


dX 


( 16 ) 


Multiplication of equation (10) by U and integration across the boundary layer together 
with the use of equation (9) result in 


dgj _ 2v oo 

dx n 3 J 0 (ay) dY 

Again multiplication and division by 6* give 

-(j<) = 

dX\ 1/ 

where 


(17) 


^*R 
U e 6* 


(18) 


j -% 


and 


2£>^ p6 /am 2 


Then by the use of equations (14) and (15) the moment of momentum equation is written as 

X2 2 

/** ^ T 

(19) 


flhL + (6*) 2 «=R 

£ dX ' ' dX 

Integration of equation (11) across the boundary layer and the use of equation (9) give 

®w (20) 


dE . 
dX U f 


V oo 

d(S/Sw) 

U e 

dY J 


iw 


Multiplication and division by 5* give 

d^T6*) 


dX 


U e 6* 


Q(Sw) 


(21) 


where 


T = 4= 
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and 


Q= 6.’ 




dY 


w 


By the use of equations (14) and (15), the energy equation can be written as 



( 22 ) 


The final equations to be solved simultaneously are equations (16), (19), and (22). 

As in reference 2 the velocity and enthalpy profiles are taken from the reverse- 
flow or "lower branch" solutions of reference 5. Because only a limited number of 
these solutions are presented in reference 5, it was necessary to obtain additional ones. 
The additional solutions were obtained by using the programed equations described in 
reference 17 for Np r = 1 and <p = 1. For these conditions the similar boundary- layer 
equations of reference 17 reduce to those of reference 5. The equations in reference 17 
were obtained by assuming similar flow and using the similarity parameters 




P w M w u e dx 


u e p dy 
^0 

f]^wVe 


for two-dimensional flow. 

The resulting boundary-layer equations are 

f' M + f(f") + /sjl + S - (f*)^] = 0 


and 


where 


S” + S’f = 0 


r x 

\ p u,„ Up dx 
_ 2 du e Jp e 

u? dx t e Pwdw 


(23) 


(24) 


(25) 

( 26 ) 
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and 

f , = u_ = JL 

u e U e 


Solutions to these equations were obtained over a range of negative values of /3 for 
t w = 0.35. As in reference 2, the velocity and enthalpy profiles are uncoupled from the 
parameter /3 by assuming that the velocity profile is related to the parameter a which 
is equal to the distance to the point of zero velocity divided by 6 ^that is, a = — 
Similarly, the enthalpy profile is assumed to be related to the parameter b which is 
equal to (that is, b = £ w ). Therefore, a reverse-flow solution to equations (25) 
and (26) for a given p and t w is assumed to be uniquely specified by the associated 
values of a and b. The integral parameters H, J, P, R, and 6* were therefore 
obtained from the reverse-flow solutions and plotted as functions of the associated values 
of a. The resulting distributions were fitted with polynomials by using the method of 
least squares. The coefficients for the polynomials are given in table I where, for 
example, 


— 2 3 4 5 

H = ci + C£a + 03a + 04a + 05a + cga 


(27) 


As the parameter E is a function of both the velocity and enthalpy profiles, it is 
a function of both a and b. Therefore, values of E were curve fitted as functions of 
a for various values of b and also curve fitted as functions of b for various values 
of a. The results of this procedure are given in tables II and HI. In table III the d 
coefficients are used to emphasize that E is a polynomial in b for the value of a 
shown. 

The parameters 6* and E shown in the tables are actually integrals over 77 
rather than Y. However, the appearance of these parameters in the final equations is 
consistent with the original definitions of Q and T, inasmuch as these parameters 
always appear as ratios to each other as in T, or 6* is multiplied by a derivative 
which can be found with respect to 77 as in Q. Because E was only available at 
discrete values of a or b, the values and derivatives of E for given values of a and 
b were obtained by linear interpolation between the available values. 

Examination of equations (16), (19), and (22) shows that the unknowns are now a, 
b, and The integral parameters are polynomials of power greater than 1, so the 

resulting set of ordinary differential equations is nonlinear. The boundary conditions 
at X = 0, which is the separation point, are = 1, a = 0, and b = (?w)x q- As this 
latter value is not generally known, it was assumed to be equal to ^ for the enthalpy 
distribution corresponding to the value of p which gives f^ = 0 (at the separation 
point). That is, the local enthalpy profile at separation was assumed to be the one asso- 
ciated with the velocity profile at separation, according to reference 5. 
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The equations for a, b, and ^6*^ as functions of X were solved numerically 
on an IBM electronic data processing machine. The result of the solution, the b dis- 
tribution, is shown in figure 8 where the ratio ^ '/t^x-O is P lotted as a function of 
X. The distribution of this ratio with x can be obtained by application of the inverse of 
the transformations. By assuming that X = 1 and using the value of (s?) obtained 
from the value of 6j at a = 0 (velocity profile parameter at separation), the heat- 
transfer rate was obtained from the expression 


q = k w 




20 



REFERENCES 


1. Carlson, Walter O.: Heat Transfer in Laminar Separated and Wake Flow Regions. 

1959 Heat Transfer Fluid Mech. Inst. (Univ. of Calif.), Stanford Univ. Press, June 
1959, pp. 140-155. 

2. Holden, M. S.: An Analytical Study of Separated Flows Induced by Shock Wave — 

Boundary Layer Interaction. Rept. No. A1-1972-A-3 (Contract NAS 5-3976), 
Cornell Aeron. Lab., Inc., Dec. 1965. 

3. Reeves, Barry L.; and Lees, Lester: Theory of Laminar Near Wake of Blunt Bodies 

in Hypersonic Flow. AIAA J., vol. 3, no. 11, Nov. 1965, pp. 2061-2074. 

4. Lees, Lester; and Reeves, Barry L.: Supersonic Separated and Reattaching Laminar 

Flows: I. General Theory and Application to Adiabatic Boundary- Layer/Shock- 
Wave Interactions. AIAA J., vol. 2, no. 11, Nov. 1964, pp. 1907-1920. 

5. Cohen, Clarence B.; and Reshotko, Eli: Similar Solutions for the Compressible 

Laminar Boundary Layer With Heat Transfer and Pressure Gradient. NACA 
Rept. 1293, 1956. (Supersedes NACA TN 3325.) 

6. Nielsen, Jack N.; Lynes, Larry L.; and Goodwin, Frederick K.: Calculation of 

Laminar Separation With Free Interaction by the Method of Integral Relations. 

Part H: Two-Dimensional Supersonic Nonadiabatic Flow and Axisymmetric 

Supersonic Adiabatic and Nonadiabatic Flows. AFFDL-TR-65-107, Pt. II, U.S. 

Air Force, Jan. 1966. (Available from DDC as AD 630 765.) 

7. Stainback, P. Calvin: Heat- Transfer Measurements at a Mach Number of 8 in the 

Vicinity of a 90° Interior Corner Alined With the Free-Stream Velocity. NASA 
TND-2417, 1964. 

8. Schaefer, William T., Jr.: Characteristics of Major Active Wind Tunnels at the 

Langley Research Center. NASA TM X-1130, 1965. 

9. Lyman, Taylor, ed.: Metals Handbook. Volume 1.- Properties and Selection of 

Metals. 8th ed., Am. Soc. Metals, c.1961, p. 423. 

10. Donaldson, I. S.: On the Separation of a Supersonic Flow at a Sharp Corner. AIAA 

J. vol. 5, no. 6, June 1967, pp. 1086-1088. 

11. Zoby, Ernest V.; and Sullivan, Edward M.: Effects of Corner Radius on Stagnation- 

Point Velocity Gradients on Blunt Axisymmetric Bodies. NASA TM X-1067, 1965. 

12. Dennard, John S.; and Spencer, Patricia B.: Ideal-Gas Tables for Oblique-Shock 

Flow Parameters in Air at Mach Numbers From 1.05 to 12.0. NASA TN D-2221, 
1964. 


21 



13. Chapman, Dean R.; Kuehn, Donald M.; and Larson, Howard K.: Investigation of 

Separated Flows in Supersonic and Subsonic Streams With Emphasis on the Effect 
of Transition. NACA Rept. 1356, 1958. (Supersedes NACA TN 3869.) 

14. Chapman, Dean R.: A Theoretical Analysis of Heat Transfer in Regions of Separated 

Flow. NACA TN 3792, 1956. 

15. Weiss, Robert F.: Base Pressure of Slender Bodies in Laminar, Hypersonic Flow. 

AIAA J., vol. 4, no. 9, Sept. 1966, pp. 1557-1559. 

16. Kruse, Robert L.: Transition and Flow Reattachment of the Separated Boundary 

Layer of an Apollo- Like Bluff Body in Free Flight at Mach Numbers 1 to 9. NASA 
TND-4645, 1968. 

17. Beckwith, Ivan E.; and Cohen, Nathaniel B.: Application of Similar Solutions to 

Calculation of Laminar Heat Transfer on Bodies With Yaw and Large Pressure 
Gradient in High-Speed Flow. NASA TN D-625, 1961. 

18. Scott, C. J.; and Eckert, E. R. G.: Heat and Mass Exchange in the Supersonic Base 

Region. Separated Flows, Pt. I, AGARD CP No. 4, May 1966, pp. 429-478. 


22 



TABLE I.- COEFFICIENTS IN POLYNOMIAL CURVE FITS 
OF INTEGRAL PARAMETERS THAT ARE 
FUNCTIONS OF a ONLY 


Parameter 


Coefficients in polynomial fit 


C 1 

c 2 

c 3 

c 4 

c 5 

c 6 

H 

0.227 

-0.283 

0.372 

-4.729 

8.420 

-2.984 

J 

.341 

-.314 

-1.051 

2.360 

-7.123 

9.357 

P 

-.036 

-.352 

-5.377 

7.716 

0 

0 

R 

1.386 

.565 

16.550 

-68.680 

212.500 

-200.300 


2.713 

1.004 

32.203 

-161.200 

481.800 

-407.200 
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TABLE n.- COEFFICIENTS IN POLYNOMIAL CURVE FITS 
OF INTEGRAL PARAMETER E(a,b) AT 
VARIOUS CONSTANT VALUES OF b 




Coefficients in polynomial fit of E(a,b) 


Value of b 

C 1 

c 2 

c 3 

c 4 

c 5 

c 6 

0.194 

-0.192 

0.431 

0.936 

3.142 

-13.810 

9.088 

.174 

-.242 

.377 

3.276 

-7.842 

9.504 

-9.533 

.164 

-.273 

.685 

- 1.000 

15.730 

-41.110 

27.860 

.155 

-.304 

.618 

.345 

10.310 

-29.940 

18.490 

.130 

-.418 

.784 

.324 

10.670 

-24.800 

8.555 

.116 

-.508 

.996 

-1.578 

20.840 

-40.330 

13.970 

.095 

-.686 

.960 

1.148 

7.246 

-.737 

-25.730 

.081 

-.846 

1.047 

2.665 

-4.069 

37.890 

-65.750 

.068 

-1.051 

1.376 

-1.716 

25.000 

-28.410 

-11.440 

.079 

-1.485 

1.200 

5.099 

-21.560 

113.600 

-147.200 

.029 

-2.177 

-.123 

31.180 

-176.400 

499.100 

-467.800 

.018 

-2.936 

.551 

22.650 

-132.800 

411.800 

-393.000 
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TABLE m.- COEFFICIENTS IN POLYNOMIAL CURVE FITS 
OF INTEGRAL PARAMETER E(a,b) AT 
VARIOUS CONSTANT VALUES OF a 


Value of a 

Coefficients in polynomial fit of E(a,b) 

d l 

d 2 

d 3 

d 4 

d 5 

d 6 

0 

-4.694 

123.00 

-1666 

12 436 

-47 483 

72 052 

.106 

-4.717 

141.30 

-2143 

17 516 

-71 831 

115 338 

.145 

-4.393 

128.20 

-1879 

14 996 

-60 586 

96 410 

.176 

-4.203 

122.20 

-1750 

13 579 

-53 457 

83 202 

.247 

-3.835 

115.70 

-1634 

12 362 

-47 576 

72 767 

.282 

-3.650 

117.20 

-1722 

13 412 

-52 817 

82 290 

.330 

-3.305 

116.50 

-1789 

14 218 

-56 364 

87 713 

.357 

-2.980 

110.50 

-1717 

13 595 

-53 396 

82 243 

.384 

-2.637 

105.80 

-1719 

14 051 

-56 634 

89 241 

.422 

-1.952 

88.22 

-1497 

12 481 

-50 788 

80 478 

.474 

-1.029 

57.51 

-1049 

9 050 

-37 471 

59 903 

.518 

-.432 

35.45 

-726 

6 711 

-29 274 

48 780 
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Face instrumentation 

Afterbody instrumentation 

Orifice number 

r 

0, deg 

Thermocouple or 
orifice number 

l ,in. (cm) 

1 

0 

1 

0.200 (0.508) 

2 

3.8 

2 

.346 (0.879) 

3 

7.6 

3 

.492 (1.250) 

4 

11.4 

4 

.638 (1.620) 

5 

15.2 

5 

.784 (1.991) 

6 

19.0 

6 

.930 (2.362) 

1 


7 

1.076 (2.733) 



8 

1.222 (3.104) 



9 

1.368 (3.475) 



10 

1.514 (3.846) 



11 

1.660 (4.216) 
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(a) Model sketches showing locations- of instrumentation. (All linear dimensions are given in inches and parenthetically in centimeters.) 


Figure L- Test configuration. 




Front view 



Oblique view 


(b) Photographs of model. L-67-8727 

Figure 1.- Concluded. 
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]-) = 4. 65 x 10 5 


Rqq -p = 6. 55 x 105 


Schlieren photographs of model flow field for various values of free-stream Reynolds number. 


L-67-8728 







(a) Typical schlieren photograph showing regions of flow field in vicinity of corner. 





Region of 
undisturbed oil 


(b) Photograph of oil flow over afterbody. L-67-8729 

Figure 3.- Photographs showing various flow regions and oil flow over afterbody. 
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Level obtained by assuming isentropic expansion to measured afterbody flow angle 


Isentropic recompression (method of ref. 13), with -S = 0.587 


Beginning of separated flow 


Rear of model 


(a) Rco D = 1.89 x iq5. 

Figure 5 .~ Afterbody pressure distributions. 
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(b) Roo d = 4 - 05 X 1[)5 - 

Figure 5.- Continued. 
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■ — Isentropic recompression (method of ref. 13), with = 0.587 
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— Level obtained by assuming isentropic expansion to measured afterbody flow angle 


Rear of model 


Beginning of separated flow 


Ic) Roon = 7.55 x 10 3 . 


Figure 5.- Concluded. 
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Figure 6.- Variation of afterbody pressure level with free-stream Reynolds number. 






Flat_plate theory 



Beginning of separated flow 


Rear of model 






(e) Rood = 4 - 05 x M 5 ; *w = 0-397; p/pj = 0.0165 for afte 


Figure 7.- Continued. 
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Flat-plate theory 

Theory, appendix A 



Beginning of separated flow 


Rear of model 



Flat-plate theory 

Theory, appendix A 



Beginning of separated flow 
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(i) Roo d = 6 * 55 * 105 ' *W = °* 39 5 


Figure 7.- 


Fiat -plate theory 

Theory, appendix A 
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p/p t = 0.015 for afterbody flow. 


Continued. 
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Theory, appendix A 
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(j) Rood = x ~ °* 386 ; P/Pt = 0.014 for afterbody flow. 
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Figure 7.- Concluded. 









